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We investigate a non-equilibrium one-dimensional model known as the raise and peel model de- 
scribing a growing surface which grows locally and has non-local desorption. For specific values of 
adsorption (u a ) and desorption (ud) rates the model shows interesting features. At u a = Ud, the 
model is described by a conformal field theory (with conformal charge c = 0) and its stationary 
probability can be mapped to the ground state of the XXZ quantum chain. Moreover, for u a > Ud, 
the model shows a phase in which the the avalanche distribution is scale invariant. In this work we 
study the surface dynamics by looking at avalanche distributions using Finite-size Scaling formalism 
and explore the effect of adding a wall to the model. The model shows the same universality for the 
cases with and without a wall for an odd number of tiles removed, but we find a new exponent in 
the presence of a wall for an even number of avalanches released. We provide new conjecture for the 
probability distribution of avalanches with a wall obtained by using exact diagonalization of small 
lattices and Monte-Carlo simulations. 
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Substrate 



FIG. 1. Three possible cases shown. (A) tile attaches to the surface, (B) tile removes a layer, (C) tile is reflected. 

I. INTRODUCTION 

The Raise and Peel Model (RPM) is a Markov process first proposed by [l| describing the evolution of a growing 
surface with a fluctuating interface in one dimension. This model has been found to belong to a new universality 
class in non-equilibrium phenomena [1 --5] . For a particular value of the adsorption (u a ) and desorption (ud) rates, the 
model exhibits a phenomenon of self-organized criticality [f| 0] where probability distributions of desorption events 
show long tails and are characterized by a varying critical exponent that depends on a single parameter given by the 
ratio of the adsorption and desorption rates [J, [5| . 

When adsorption and desorption rates are equal, the model becomes solvable. This goes back to a connection 
established by Razumov and Stroganov pj-fl0| which relates the two-dimensional dense 0(n = 1) fully packed loop 
models (enumerating the stationary state probability distributions of RPM) to those a ground state wavefunctions of 
the XXZ chain with L sites [3, HH, HIl- Moreover, the spectra can be obtained by conformal field theory with charge 
c=0 [13, Hil- This offers a nice mathematical structure, which allows to make conjectures using small lattices for 
expressions of physical quantities that remain valid for quantities for any system size. 

In this work, we study the effect on avalanches in the presence of a wall (RPMW) since little is known about this 
effect when the boundary is allowed to fluctuate. Some other interesting results with the wall have been reported 
for example in [ll| ■ In Section [II] we describe the stochastic rules for the model with and without a wall and 
highlight some of the known results for these two cases. In Section Hill we compare the energy spectra of the stochastic 
Hamiltonian of the XXZ quantum chain for different spin sectors. Lastly in Section ITVl we compute critical exponents 
for avalanche distributions for RPM and RPMW and derive new conjectures for probability expression with a wall. 

II. RAISE AND PEEL MODELS 

The Raise and Peel Model (RPM) describes a growing and fluctuating interface. An initial configuration is chosen 
and tiles are dropped onto the surface with a certain probability. Three different processes can happen as shown in 
figure [TJ With some probability Pi = L+ 1 a _ 1 a tile lands in site i = 1, L — 1, while in RPMW the site i=0 is chosen 
with some probability Pq = ■ Depending on the slope of surface at the i th site, one of three things can occur. 

• Case A: Tile hits a local minimum 

if site i > 0, with some probability u a the tile attaches to the substrate, else if site i=0 is chosen then half-a-tile 
attaches to the boundary with rate 1. 

• Case B: Tile hits a slope 

With Probability Uf, the tile peels off tiles within a cluster such as the local local height at every site in the 
cluster decreases by two: hi —¥ hi — 2 (a tile has a height of 2); In other words, a tile may only remove one layer 
of tiles above the point of contact. 

• Case C: Tile hits a local maxima 
Tile reflects and nothing happens. 

This stochastic process in continuum time is given by the master equation [ItI Il8j 

±P a {t) = -Y J H a ^Pp{t) (1) 



where P a (t) is the (unnormalized) probability of finding the system in one of the states |a) at time t, and H a p is 
the rate for the transition |ev) — > \j3). Since this is an intensity matrix, there is at least one zero eigenvalue Q and its 
corresponding eigenvector |0) gives the probabilities in the stationary state 
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(o|fr = o, (o| = (i,i,...,i) (2) 

H\Q) = 0, |0) = VP a |a), P a = lim P a (t) (3) 

* — 4 £— >-oc 
a 

For the special case in which the rates are equal u a = iid, the Hamiltonians can be written in terms of the Temperley- 
Lieb algebra defined in terms of generators ei, . . . , satisfying the following commutation relations fl9l | 

e i = e i 

[e i ,e j ]=Q for \i-j\>2 (4) 



while the one-boundary term eo at site i = satisfies the following constraint |20t l21j 

e" = e , eie ei=ei, e ej = eje if i > 1 (5) 

There are many representation to the Temperley-Lieb algebra including "blob-algebra" . For our purposes it becomes 
convenient to view the generators in terms of a tile and half-tile at the boundary. 








Products of generators at different sites, following the algebra relations, reduce the products to a subset of unique 
configurations. This is illustrated by the following two examples. 



e4e 6 eie3e 5 e7 = 




012345678 




In terms of these generators, the stochastic Hamiltonians H for RPM and for RPMW with a rate a at the 
boundary can be written as: 



L-l L-l 



H = ^2(l-ei) = a (l-e ) + ^(l-O (9) 



i=l 



The ground state eigenvectors of the intensity matrices (Eqn. [5]) have remarkable combinatorial properties [3, d, [l(| ■ 
The normalization on the ground state eigenvectors will be used in section [V] to derive conjectures for probabilities. 
The following two examples illustrate the combinatorial properties for a small lattice with L=6 in RPM and for L=4 
in RPMW. 



H = 



/ -2 2 2 2 \ 
'1-30 1 * 
1 0-3 1 
1 1-3 2 
V 1 -4 / 



#(<*) = 



/ -3 a 

I- 2 1 a 

II- 21 
1 1 -2 -a 
1 

\ 



\ 




a 


1 -1-a 2 
1 -2 - a / 



The wave functions |0) (|0)^ a ^) normalized to have the smallest entry equal to 1 (or a) are given by: 
|0) = (11,5,5,4,1) |0) (a) = (a 2 ,3a(2 + a),2a(3 + a),3a,3(2 + a),3) 

The normalization factor (0|0) for RPM is then given by (L=2n) 

(0|0) = A v (2n + 1) = 1, 3, 26, 646, ... 
while the normalization factors (0|0) for RPMW with a=l is given by 

(1) (0|0) (1) = N s (2n)A v {2n+ 1) = 2,33,26,4420,... 
where Av(2n + 1) is the number of vertically symmetric (2n + 1) x (2n + 1) alternating sign matrices [221424 

A v {2n + 1) = Dy^(3j + 2) ^±l»±|j = 1, 3, 26, 646, ■ ■ • 

and Ns(2n) is the number of cyclically symmetric transpose complement plane partitions [2~ll [22j | 



7v 8 (2n) = iK^(3j + i) ^, = 1, 2, n, 170, 7429, • • 



(4j)!(4j + 1)! 



(10) 



(11) 



(12) 



(13) 



(14) 



III. ENERGY SPECTRA AND SPACE-TIME PHENOMENA 



For u = Ua/ud = 1, the finite-size corrections to the ener gy spectra of the intensity matrices H and flw are given 
by a conformal field theory with central charge (c = 0) [25H27J 

E n =E + 7rv{As + k) +0(L- 1 ) keZa,2dv = ^- (15) 



A s = £^ = 0,0,i,l... S = 0,i,l,|,... (16) 

In table U the excited energy states E n for different spin sectors are compared to the numerical estimations obtained 
by diagonalizing the intensity matrices H and for RPM and RPMW (Eqn. [9]). 



TABLE I. Excited energy states in units of ^ for different spin sectors, and numerical approximations to E n using L=18 for 
RPM and L=16 RPMW 



A s + k 


1 2 3 4 5 


s = 

« = 5 
a = l©« = l/2 


2 3 4 4 5 
1 2 3 3 4 
1 2 2 3 3 


n 


1 2 3 4 5 


RPM 
RPMW 


0.0000 2.0009 3.0035 4.0247 4.0257 5.037 
0.0000 1.0015 2.0087 1.9954 3.0030 3.0158 
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The functional dependence of the time evolution of distributions can be predicted. The expectation value of 
observables can be described by stochastic dynamics as: 

(X)(t) = (0\Xe- m \*(0)) (17) 

where the initial state ^(O)) can be expanded in a complete eigenbasis characterizing the system: |\&(0)) = J2 n c n\^n), 
and H is the stochastic matrix describing the system. Since H is an intensity matrix, the lowest eigenvalue is zero, 
hence the lowest non-zero eigenvalue E\ is expected to dominate the time evolution for large times. 

Let us consider the effect of adding a wall on temporal profiles of quantities describing the system. The temporal 
average of a quantity will be denoted by the Family- Vicsek [28| scaling form: 

X <M-1$- 1 ~ X & ,18 » 

The average number of clusters K(x) = &hi,o) is plotted in Fig. [2] in the form given by Eqn. fTS)) . The data 
collapse shows that the average number of clusters K(t,L) has a critical exponent given by z — 1 for different values 
of the rate a: (a=0 RPM, a=l RPMW). The long time decay is described well by the exponential given by E\ = 
for RPM and by E 1 = ^ for RPMW as is expected from Eqn. (jT7|) for large times. 




FIG. 2. The time evolution of the average number of clusters K(t) for RPM and RPMW. The lines are the expected decays 
with K{t, L) oc e~ Elt/L with £i = ^ (RPM) and Ei = ^ (RPMW) 



IV. AVALANCHES 





012345678 012345678 
FIG. 3. Example of an even and odd avalanche in RPMW. Avalanches occurring on RPM release only an odd number of tiles 



The raise and peel model exhibits events where layers are evaporated from the substrate when a tile from the gas 
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TABLE II. Estimates for the critical exponents in S(v, L) 
and L'=8192 



r F(v/L D ) for even and odd avalanches using lattices L=4096 



D 



1/u 



1.0 



0.45 



0.005 



RPM [4] odd 1.004 1.026 1.006 

RPM[this work] odd 0.992 ± 0.058 1.015 ± 0.006 1.006 ± 0.0001 

RPMW[this work] odd 0.994 ±0.025 1.008 ±0.002 1.006 ±0.0001 

RPMW[this work] even 0.980 ±0.117 1.013 ±0.008 1.006 ±0.0001 



RPM [4] odd 3.000 2.25 2.00 

RPM[this work] odd 2.977 ± 0.079 2.224 ± 0.016 2.011 ± 0.001 
RPMW[this work] odd 3.003 ±0.071 2.237 ±0.012 2.011 ±0.001 
RPMW[this work] even 1.932 ± 0.222 1.280 ± 0.074 1.020 ±0.011 



hits the interface. The number of tiles removed defines the size of an avalanche. While this number is always an odd 
number in RPM, in RPWM there is the possibility of an even number of tiles removed whenever an avalanche touches 
the boundary. This is illustrated in Figure [3] 

It is known that the raise and peel model exhibits self-organized criticality @, 0] in the regime for u > 1 [![■ 
Desorption processes being non-local results in avalanches lacking a characteristic length-scale. Their distribution 
S(v,L) therefore appears as a power-law which might be described in the finite-size scaling (FSS) form 0, H 

S(v,L)=v- r F(^ 5 ) (19) 
In order to obtain the exponents, the method of moments is used [1, [29[. Using the scaling form (Eqn. [T9|) we have: 

(v m ) L = { S{v,L)v m dv (20) 



= J v- T F{j^)v m dv (21) 
= J W - T L- DT F(w)w m w mD L D dw (22) 
= £f (l+m-T) f w m -TF{w)dw (23) 



r m (24) 

where we have used w = v/L D to get the scaling dependence with L. We can get an estimate for the exponent by 
looking at the ratio: 

(v m ) L /{v m ) L , = (L/L>y (m) (25) 
and in this manner the exponent a(m) can be estimated as [l|: 



= ln((v m ) L /{v m ) L .) = f for m < r - 1 

1 ' ln(L/L') [ D(l + m-r) for m>r-l 1 ' 

A linear fit to Eqn. ([2"6"f for m > r — 1 gives an estimate for the values of D and r. To get an idea of the spread of 
these values we ran several Monte-Carlo simulations to find the variation of the distribution resulting from different 
seeds. The results are shown in table [TT1 

Results in this table show that the critical exponents for an odd number of tiles removed remains unchanged by 
adding a wall. However, we found that for an even number of tiles the power-law exponent r decreases by about one. 
We also found an interesting effect on the finite-size scaling function with the addition of the wall. Fig. 0] shows the 
scaling function for RPM and RPMW for u=l, where we use r = 3.0 and r = 2.0 respectively, for odd number of 
tiles and even number of tiles, while D is kept fixed at D = 1. Fig. [5] shows a similar plot for l/u=0.005. The data 
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collapse for large lattices on these plots confirms the FSS form (Eqn. H9)) . 

V. CONJECTURES FOR PROBABILITIES 

Simple conjectures for the probabilities of absorption, desorption, and reflection can be written down by considering 
the rate of change between the different states. The probability to loose (or gain) v tiles can be written as [l|: 

P(v,L) = J2 5(v(v')-v(v)-y)w, r ^P, r /(0\0) (27) 

where P l; //(0|0) is normalized probability to be in the state \v(rj')} (see Eqn. 1 10[) and Wrf-^ is the transition rate 
w-q'^n from state \v(r)')) to state \v(i~i)). For u=l , the normalization (0|0) is given by Eqn. ([TT|) and Eqn. (fl~2"l) for RPM 
and RPMW, respectively. The probability of absorption (P a ) for a given L has been given by Alcaraz et al for RPMW 
[l5j |. In this section we will present new conjectures for the probabilities of desorption (Pa) and reflection (P r ) for 
RPMW. Probabilities for RPM where first reported in [lj. These expressions results in quotients of parabolas as seen 
in table Hill Notice the denominator for the probabilities are different with the addition of the wall since normalization 
expressions for (0|0) (Eqns. [TT1 and [1~2")) are different for the cases with and without a wall [l5j . 

TABLE III. Conjectures for the probabilities of an absorption event P a = P( — 1,L), desorption event Pd = P(v > 0, L), 
reflection event P r = P(0, L) 



Prob.[l, 15] 


Pa 


Pd 


Pr 


RPM (v G odd ) 


3L(L-2) 


2(L-2)(L+2) 




4(2L + 1)(L-1) 


4(2L+1)(L-1) 




RPM (v G even ) 






3L 2 +2L+4 
4(2L+1)(L-1) 


RPMW (v G odd ) 


6L 2 +8L-5 


4L 2 +5L+9 




4(2L + l)(2L+3) 


4(2L+l)(2L+3) 




RPMW (v G even) 




1.795L-4.562 


6L' 2 + 17.211L+12.50 




4(2L+l)(2L+3) 


4(2L+l)(2L+3) 



It is interesting to note that although this is a non-equilibrium system, simple expressions for probabilities can be 
obtained. The mean size of an avalanche (v) in the stationary state can be conjectured to be given by a mean-field 
expression 




v/L v/L 

FIG. 4. The scaling function F(~) with rates l/u = 1.0 for v £ odd number with r = 3.0 and D — 1.0 (LEFT) and v £ even 
number with r = 2.0 and D — 1.0 
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L = 2048 
L = 1024 



1/u = 0.005 
RPM 




0.(14 




L 


= 2048 






■■ L 


= 1024 




0.2 0.4 



v/L 



0.6 0.8 



FIG. 5. The scaling function F(£) with rates 1/u = 0.005 for v G odd number with r = 2.0 and D = 1.0 (LEFT) and v € even 
number with r = 1.0 and D = 1.0 



y) L = P a (L)/P d (L) 



(28) 



A 
> 

V 



1.55 



1.45 



1.35 



1.25 



<v>(L) = v_(1 + B/L + 




RPM 
RPMW 

s — 



j i i i 



<v>(L)/v - 1 )*L 



- 1 — i — i i i i 1 1 1 1 — i — i i i i i | 1 — 

j -B-B -B-B EEB5E3- - - S" ■ 



T 1 — I I I I I 



v o cocooro — e — D 

_l i I i i i I i i jjj iii 



100 



1000 



10000 



—I l_ 



_l I I I I 



_l I l_ 



1000 



FIG. 6. Average avalanche size (v)(L) plotted from the ratio of P a (L)/Pd(L) (lines) compared to Monte Carlo data for RPM 
(O) and RPMW (□) 
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This is a quotient of parabolas, and in the large limit (L 3> 1) 

Pa{L) _ aL 2 + BL + 1 
P d {L) aL 2 + bL + c 

aL 2 + /3L + 7 



a 



(2 + _^+0(L- 2 ))(l-±+0(L- 2 )) 
a aL aL 

,ol 8 ab ,„, o, s 

- + ^-^ + or 2 

a a L 

,a 



F + >--)) 

a aL a 
a _L a a 



(29) 



We assume that a ^ and a ^ and dropped terms of the order of 0(L 2 )). From Table Hill we see that for RPM 
we have 



P d (L) 2' L v 3 2' 

(30) 



while for RPMW we have: 



P a (L) 6 1 8 6.795 
P d (i) 4 l L l 6 4 JJ 
3, 0.365, 



(31) 



The estimated average from the ratio of the conjectures agrees quite well as can be seen in Fig[6l This remarkable 
since this describes a non-equilibrium system. 

It is interesting to see that the leading term on this expansion is universal, whereas the correction term depends on 
the details of the model (i.e. whether there is a wall or not). This is a similar behavior as in finite-size scaling of the 
concentration of particles in certain reaction-diffusion systems (30l - [32l | and surface exponent corrections to quantum 
chains using different types of boundary conditions [33l435| . 

The simple functional form for the probabilities shown in table IIIII suggests we can guess a general quadratic 
expression in L for the probabilities P(v,L) by fixing the denominator as: 

p / n a(v)L 2 + b[v)L + c{v) 

PRPM{V > L) = 4(2L + !)(£-!) (32) 



p / n a(v)L 2 + b(v)L + c{v) 

Phpmw(v, L) = 4(2£ + 3)(2L + 1) (33) 

and fitting for the parameters {a(v), b(v), c(v)} in the forms (|32|) and ([33]) . The behavior of these parameters with 
respect to v is shown in Fig [71 As expected, the quadratic term shows a power-law behavior cx v~ 3 . The linear 
and constant terms however do show different behavior, but we were not able to reduce it to an analytical form. 
Consistency on the fits demands that (see table 



(34) 
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FIG. 7. Behavior of the quadratic term a(v), linear term b(v) and constant term c(v) obtained from fits to Egn. l32l and Egn. 1331 



^Phpmw(v,L) = 4(2L + 3)(2L + 1) (35) 

The sums over the parameters J2 V { a i v )>b(v), c(v)} are shown in Table ITVl for RPM and RPMW. We see that the 
quadratic a(v) and linear b(v) terms capture the behavior in v quite well, whereas as shown in Fig [7] the constant 
term c(v) is dominated by fluctuations on the fits. 





£„&(«) J2 v c(v) 


RPM 
RPMW 


1.99999 -0.00375 -13.8002 
3.99997 4.96604 30.4445 



TABLE IV. Sums over v for the parameters resulting from the fits 



VI. CONCLUSION 



We studied the non-equilibrium statistical model known as the raise and peel model. We have confirmed that 
this model retains several features as predicted from conformal invariance for stochastic profiles characterizing the 
system when changing the boundary conditions. We allowed one boundary to fluctuate and demonstrated that the 
temporal profile of a stochastic quantity follows its expected behavior from stochastic dynamics where the long time 
behavior is dominated by the lowest non-zero eigenvalue. We study the surface dynamics by looking at avalanche 
distributions exhibiting power-law distributions. Using the finite-size scaling formalism we confirm the universality 
exponent r = 3.0 for the Raise and Peel with different boundary conditions and we identified an even/odd effect with 
a new exponent r = 2.0 for avalanches with an even number of tiles removed. We also found new conjectures for 
the probability of desorption and reflection with a wall added to the system and checked that they agree with Monte 
Carlo data. 
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